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ABSTRACT 



o 

r-j I Context. Recent work based on a global measurement of the ICM properties finds evidence for an increase in the iron abundance in 

galaxy clusters with temperatures around 2-4 keV up to a value about 3 times higher than is typical of very hot clusters Zpe — 0.25Zq . 
I . Aims. We have started a study of the metal distribution in nearby X~ray luminous poor galaxy clusters in the temperature range 3-4 

O ' keV with measured high abundances. Our study aims at spatially resolving the metal content of the central regions of the ICM, in 

particular characterizing how our measurements are biased by the diagnostics adopted and by the possible multi-temperature nature 
C/5 , of the projected observed spectra. We report here on a 42ks XMM-Newton observation of the first object in the sample, the cluster 

C^ ,. Abell 2028. 

Methods. We selected interesting regions of the ICM to analyze the spatially resolved structure of projected temperature and abun- 
dance, thereby producing two-dimensional maps. 

Results. The XMM-Newton observation of the first object of the sample, the cluster Abell 2028, reveals the complex structure of the 
cluster over a scale of ~ 300 kpc, showing an interaction between two subclusters in cometary-like configurations. Cold fronts have 
been detected at the leading edges of the two substructures. The core of the main subcluster is likely hosting a cool corona. We show 
that a one-component fit for this region returns a biased high metallicity. The inverse iron bias is caused by the behavior of the fitting 
code in shaping the Fe-L complex. In the presence of a multi-temperature structure of the ICM, the best-fit metallicity turns out to be 
artificially higher when the projected spectrum is modeled with a single temperature component, while it is not related to the presence 
' of both Fe-L and Fe-K emission lines in the spectrum. After accounting for the inverse iron bias, the overall abundance of the cluster 

' is consistent with the one typical of hotter, more massive clusters. 

, Conclusions. We caution against interpreting high abundances inferred when fitting a single thermal component to spectra derived 

from relatively large apertures in 3-4 keV clusters, because the inverse iron bias can be present. Most of the inferences trying to relate 
high abundances in 3-4 keV clusters to fundamental physical processes will probably have to be revised. 



Key words. X-rays: galaxies: cluster: general - intergalactic medium - galaxies: clusters: individual: Abell 2028 - galax- 
ies:abundances 

1. Introduction sequence of the past history of star formation in the cluster 

galaxies and of the processes responsible for injecting enriched 
Clusters of galaxies have established themselves as a unique en- ^^^^ (^e ICM. (For recent review s of the o bservation al and 

vironment for measurmg elemental abundances and for study- theoretical aspects. seelWerner et al.i r2008: Bor gani et al.ll2058l: 



ing the chemical enrichment history of the Universe. Their large ISchindler & Diaferioll2008l andT^ences therein.) 
potential well retains all the metals produced by their mem- 
ber galaxies, allowing in principle a unified description of the Most of the appeal of the X-rays determining elemental 
thermodynamical properties of the hot phase (the diffuse X- abundances in the hot gas stems from the apparent robustness of 
ray emitting gas) and of the cold phase (the cluster galaxies), the measurement: most of the observed emission lines arise from 
Indeed, the distribution of metals in the ICM is a direct con- well understood hydrogen- and helium- like ions and their equiv- 
alent widths can be, under the reasonable assumption of colli- 
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the reliability of the Fe-L shell modeling (e.g. iRenzin il [19971 
lArimoto et al .''l997). it has been shown that this is not a par- 
ticular issue (Buote & Fabian 1998; Buote et al. 2003b) 

Unfortunately this intrinsic simplicity of the measurement 
has faced the limitations of the X-ray satellites in terms of 
the shape of the instrumental response, bandpass, and spectral 
and spatial resolution leading to observational biases. Correct 
modeling of the temperature structure is crucial, in particu- 
lar when dealing with spectra extracted from a large aperture 
centered on the core of clusters and groups. Given the pres- 
ence of strong and opposite gradients in the temperature and 
metallicity profiles (i.e., coo ler regions are more metal rich, see 
iDe Grandi & Molendil l200ll) . in the spectra of central regions 
photons coming from regions of different temperature and abun- 
dances are mixed together. These objects and these regions natu- 
rally draw the attention of X-ray observers because of their high 
surface brightness, thus delivering high signal-to-noise (S/N) 
spectra. 

The first important recognition of a bias in the measure- 
ment of elemental abundances was the description of the "iron 
bias" dluote & Fabian 1998; Buote 2000a b). Many ASCA 
studies found significant subsolar values f or the iron abun- 
dance in grou ps and elliptical galaxies (e.g. iDavis et anil999t 
iMatsumoto e t al. 1997), generally less than the stellar values in 
these systems and less than those observed in galaxy clusters. 
iBuotd J2000b) demonstrated that this was a bias resulting from 
fitting with a simple single-temperature a multi-temperature 
plasma, resulting in best-fitting elemental abundances that are 
too low. The typical ASCA spectrum for these objects was ex- 
tracted from a large aperture in the core because of the large 
point spread function (PSF) of its telescope. Since the descrip- 
tion of this bias, and with the help of the quantum leap in- 
crease in the quality of data with XMM-Newton and Chandra, 
this kind of error has been recognized and made less severe. The 
best-fitted abundances in the multi-phase cooling core regions 
of clusters and in groups have been shown to increase when 
two-temperat ure models or more complicated DEM models 
are used (e.g. Molendi & Gastaldello"200l'; 'Buote et al."2003b^, 
Humchrev & Buote 2006; Werner et al. 2006; Matsushita et al. 
20071) . These new determinations helped to mitigate some of the 



trend seen in a plot of iron-mass-to-light ratios (IMLR) versus 
cluster temperature, which seem to drop by almost three orders 
of magnitude below 1 keV (iRenzinill 19971) . 

Another trend in the ICM abundance versus cluster tem- 
perature has been shown with increasing evidence in recent 
years. Using the ASCA archive observations of 273 objects 
iBaumgartner et al] (l2005h showed that clusters with gas temper- 
ature between 2 and 4 keV have a typical mean abundance that 
is larger by up to a factor of 3 than hotter systems . Investigating 
a sam ple of 56 clusters with Chandra at z > 0.3. lBalestra et alJ 
( 1200 7) confirms the trend of the Fe abundance with cluster tem- 
perature also in the higher redshift clusters. The observed trend 
may reflect a more efficient star formation in smaller clusters, 
as also suggeste d by optical and near-infrared observ ations of 
nearby systems (iLin et all 120031: iGonzalez et al.ll2007h . On the 
contrary it may be argued that it is difficult to conceive a phys- 
ical mechanism responsible for making clusters in the 3-4 keV 
range so unique a mas s scale for chemical enrichment (see for 
example discussions in iRenzinll 1 997l l2004l) and that some plau- 
sible measurement bias might be in action, which has to work in 
the opposite way to the know Fe bias; i.e, it has to bias iron abun- 
da nces hi^h. A funda mental step has been made in this direction 
bv lRasia et al.l ( l2008l) . who analyzed mock XMM-Newton obser- 
vations of simulated galaxy clusters finding a systematic over- 



estimate of iron for systems in the 2-3 keV range. The reason 
they put forward for this discrepancy is the complex tempera- 
ture structure of the simulated clusters, due to projection and 
low spatial resolu t ion, an d to their particular temperature range. 
ISimionescu et al.l (l2009l) support this explanation by analyzing 
a deep XMM-Newton observation of the high-luminosity cluster 
Hydra A finding a biased high Fe abundance at a level of 35% in 
the central three arcminutes, for the first time dubbing this kind 
of bias as the "inverse" Fe bias. This appealing solution has been 
suggested by'Rasia et aL (l2008h as a possible explanation of the 
Baumgartner et al. (20051) results. 

With the aim of shedding more light on these issues by go- 
ing beyond single-aperture measurements, we started with the 
present work a study of the metal distribution in iron-rich poor 
galaxy clusters ( 3 < fcTgas < 5 k eV) at redshift below 0.1. 
From the work of lBaumgartner etaP (f2005), we selected the ob- 
jects whose measurements lie at the highest end of the metal 
distribution. Here we present the analysis of the XMM-Newton 
observation of Abell 2028, the first object in the sample. 

Abell 2028 is a Bautz-Morgan type II-III galaxy cluster of 
Richness Class 1 located at redshift 0.0777 ( Struble & Rool 
1999). At this redshift 1 arcmin corresponds to 88 kpc. The 
Galactic absorption in the cluster's direction is 2. 46 x 10^° cm^^ 
dKalberla et alj|2005h . iBaumgartner et alj (l2005h used ASCA to 
measure within an aperture of 7' an emission-weighted temper- 
ature of 3.91 ± 0.61 keV and metaUicity of 0.77t J'^^Z© (error 
at the 90% c.l.). 

The outline of our work is the following. In Section 2, we 
describe the observation and data reduction. In Section 3 we 
describe the morphology of the cluster and discuss the surface 
brightness profiles in selected regions. In Section 4 we discuss 
the spectral analysis and the temperature and abundances maps, 
highlighting the relevance of the inverse iron bias in the cen- 
tral region of A 2028. We discuss the results of our analysis in 
Section 5. We summarize our results and draw conclusions for 
the present study in Section 6. We show the results of simulations 
aiming at characterizing the inverse iron bias in the Appendix. 
Throughout this work, if not otherwise stated, we plot and tabu- 
late values with errors quoted at the 68.3 per cent (Icr) level of 
confidence, assuming a Hubble constant = 70hjQ km s~^ 
Mpc^^ and flm = 1 — ~ 0.3. All the metallicity estimates 
refer to the solar abundance in Anders & Grevesse (1989) just 
for ease of comparison with previous work, because this set of 
solar abundances is now outdated. 



2. Observation and data reduction 

The cluster Abell 2028 was observed by XMM-Newton on July 
29, 2007. The data were reduced with SAS v8.0 using the tasks 
emchain and epchain. We only considered event patterns 0-12 
for MOS and for pn, and the data were cleaned using the stan- 
dard procedures for bright pixels and hot columns removal (by 
applying the expression FLAG == 0) and pn out-of-time cor- 
rection. The energy scale of the pn over the whole spectral band- 
pass was further improved with the task epreject. Periods of high 
backgrounds due to soft protons were filtered using a threshold 
cut method by means of a Gaussian fit to the peak of the his- 
togram of the time b i ns of the hght curve (see Appendix A of 
iPratt & Arnaudll2002l:lDe"Luca & Molendil2004 ) and excludede 
where the count rate lies more than 3a away from the mean. The 
light curves were extracted from regions of least source contam- 
ination (excising the bright object core in the central 8' and the 
point source list from the SOC pipeline, after visual inspection) 
in two different energy bands: a hard band, 10-12 keV for MOS 
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Fig. 1. Mosaic of the MOSl, MOS2, and pn images in the 0.5-2 
keV energy band smoothed on a 12" scale. The image has been 
divided by the summed exposure maps to correct for exposure 
variations. 



and 10-13 keV forpn (using 100s bins), and a wider band, 0.5-10 
keV (using 1 0s bins), as a safety check for possible flares with 
soft spectra dNevalainen et al.ll2005l: iPradas & Kerpll2005h . The 
flaring periods thus determined were further checked by visual 
inspection of the light curves, resulting in net exposures of 40.2 
ks, 38.9 ks, and 29.5 for respectively the MOSl, MOS2, and pn 
detectors. Point sources were detected using the task ewavelet 
in the energy band 0.5-10 keV and checked by eye on images 
generated for each detector. Detected point sources from all de- 
tectors were merged, and the events in the corresponding regions 
were removed from the event list, using circular regions of 25" 
radius centered at the source position. The area lost by point 
source exclusion, CCD gaps, and bad pixels was calculated using 
a mask image. Redistribution matrix files (RMFs) and ancillary 
response files (ARFs) were generated with the SAS tasks rmf- 
gen and arfgen, the latter in extended source mode. Appropriate 
flux-weighting was performed for RMFs, using our own dedi- 
cated software, and for ARFs using exposure-corrected images 
of the source as detector maps (with pixel size of 1 ', the mini- 
mum scale modeled by arfgeri) to s ample the variation i n emis - 
sion, following the prescription of ISaxton & Siddiquil (l2002h . 
For the background subtraction, we adopted a complete model- 
ing of the various background components as described in de- 
tail in iGastaldell o et alJ (120071) . We accurately model the ac- 
tual sky background using a region free of cluster emission at 
r > 10' to the east where the signal-to-noise ratio of the sur- 
face brightness profile in the 0.5-2 keV band falls below 2. Abell 
2028 is projected towards a region of enhanced Galactic back- 
ground, and in agreement with previous studies (see for exam- 
ple the case of Abell 2029, a cluster in the same region of the 
sky in IVikhlinin et al.l l2005h . we found it necessary to add a 
third thermal component at ~ 0.4 keV. We checked the obser- 
vation for contamination by solar wind charge exchange. The 
ACE SWEPAIVQ proton flux was less than 4 x 10^ protons s ^ 
cm~^ which is typical of th e quiescent Sun (e.g.JSnowden etaP 
l2004t iFuiimotoet all 120071: lYoshino etal]l2009h . and die fight 
curve accumulated in a 0.5-2 keV energy band (after soft proton 
flares removal) did not show any variability. The spectra of the 



out-of-field-of-view events of CCD 4 and 5 of MOSl and CCD 
5 of MOS 2 showed an anomalously high flux in the soft band 
(see Kuntz & Sno wden,2008.) so they were therefore excluded 
from our analysis. 



3. X-ray image and surface brightness profiles 

For each detector, we created images in the 0.5-2 keV band and 
exposure maps at the monochromatic energy of 1 keV and then 
combined the images into a single exposure-corrected image, 
smoothed on a scale of 12", as shown in Fig[T] The X-ray emis- 
sion clearly reveals two merging subclusters, a more luminous 
SB component (which we will call subcluster A in the follow- 
ing and labeled accordingly in Fig|2]i and a less luminous NW 
one (subcluster B in the rest of the paper and labeled accord- 
ingly in Fig|2l). The "cometary" shape of the two subclusters 
seems to suggest the direction of the merging, with the sub- 
cluster A coming from the NE direction and the subcluster B 
coming from SW. A further confirmation of this interpretation 
comes from the presence of surface brightness jumps in the sug- 
gested direction of motion of the two subclumps: we extracted 
from the combined MOS and pn exposure-corrected images pro- 
files from selected angular sectors as depicted in Fig |2] centered 
on the surface brightness peaks of the two subclusters, with bins 
accumulated to have 100 counts. The surface brightness jumps 
are visible in the SW direction for subcluster A (sector Al, cen- 
ter panel of FiglJI in the radial range 40"-70" and in the NE 
direction for the subcluster B (sector B 1 , right panel of Fig|2]i in 
the radial range 30"-70". When coupled with the spectral infor- 
mation (see Section lOT i. these surface brightness jumps can be 
interpreted as merging cold fronts. For comparison plots of the 
profile in the direction of the tails of the two subclusters are also 
shown (sector A2 for subcluster A and sector B2 for subcluster 
B), highlighting the clear asymmetry in the surface brightness 
distribution of the two subclusters. 

The X-ray peak location of the main subcluster A is consis- 
tent with the position of its brightest central galaxy, GIN 416. In 
Fig|3]we show the surface brightness profiles extracted from the 
combined MOS image and accumulated in bins of 1.1" (corre- 
sponding to the MOS physical pixel size; MOS is better suited 
to this type of analysis because it oversamples the core of the 
PSF, on the contrary the pn pixel physical size is 4.1"). In the 
same figure we plot the profile of the bright point source (RA 
15:09:34, DEC H-07:30:54) SE of the main subcluster and the 
XMM-Newton profile of Abefl 2029 (OBSID 01 1 1270201), one 
of the clusters with the most peaked surface brightness profile 
and at the same redshift of A2028. Profiles are scaled to match 
their peaks of surface brightness. In the innermost 8", the sur- 
face brightness data indicate of an excess. We fitted the inner 
30" of the A 2028 A profile with a Gaussian plus a power law, 
the latter component modeling the ICM emission. The best fit 
model, with a x^/dof = 29/28 gives a width for the Gaussian 
of o- = (4.9 ± 0.3)", which corresponds to an FWHM of 
(11.5 ± 0.7)" bigger than the on-axis FWHM of the XMM- 
Newton PSF, which is 4.4" for MOS at 1.5 keV. This excess 
central component, which is more extended than a point source, 
may be a cor ona associated with the BCG, as seen in m any merg- 
ing clusters (IVikhlinin et all l2001al: fSun et alj 1200'^) . Coronae 
are sources extended on kpc scales, and they are the remnants 
of galactic X-ray coronae with temperatures around 1 keV, com- 



the data available at|http://www.srl.caltech.edu/ACE/ACS/| 



" XMM-Newton Users Handbook, 

[http://xnun.esac.esa.mt/external/xrnin_user_support/docurn entation/uhb/node 1 7 
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Fig. 2. Left panel: Exposure-corrected MOS image in the 0.5-2 keV band with the angular sectors used to extract the surface brightness profiles 
shown in the other panels. Point sources were masked when extracting the profiles. The positions of the surface brightness jumps discussed in the 
text and shown in the next panel are indicated by arrows. Center panel: Surface brightness profiles for the selected angular sectors shown in the 
left panel for the main subcluster A. The radial ranges where surface brightness profiles change rapidly because of the surface brightness jumps 
are indicated by the shaded gray area. Right panel: same as the center panel, but for the smaller subcluster B. 



4. Spectral analysis 
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Fig. 3. Surface brightness profile of the core of subcluster A of 
A2028 compared to the one of the bright point source SE of 
the subcluster and to the A2029 cluster, one of the clusters with 
the most peaked surface brightness profile. The profiles have 
been scaled to match the peak of the radial profile of A2028 A: 
the point source profile was multiplied by 0.29 and the A2029 
profile by 0.058. The dotted red line corresponds to the best-fit 
GaussianH-power-law model of the profile of the core of A 2028 
A discussed in the text. 



pressed by the hot gas of the cluster The spectral information 
we present in Section l4!2l will strengthen this interpretation. 



4. 1. Temperature and abundance maps 



We extracted spectra in a series of regions as depicted in the 
left panel of FigH) chosen to select relevant regions of the ICM. 
Given the complex geometry and the relatively low number of 
counts we did not de-project the spectra and the results quoted 
in this paper come from fitting projected spectra. They consist in 
two concentric annular regions centered on the emission peak of 
the subcluster A (RA 15:09:28 DEC H-07:33:20) with bounding 
radii 0'-0.5', 0.5'- 1' and four sectors with bounding radii l'-2' 
and 2'-4'; a circle of 0.5' radius centered on the emission peak 
of the subcluster B (RA 15:09:20 DEC H-07:38:15), two sectors 
with bounding radii 0.5'-l' and a box of 4'x2.4' encompass- 
ing the tail. The radii of regions 2 and 3 and of regions 7 and 8 
have been determined to match the surface brightness jumps dis- 
cussed in Section[3]as closely as possible, in order to investigate 
the temperature variation across these fronts. We fitted the back- 
ground subtracted spectra with a single APEC thermal plasma 
model (Smith et al. 2001) in the 0.5-10 keV band with the ab- 
sorbing column density fixed at the Galactic value. The tem- 
perature and abundance values of the regions of the map thus 
obtained are shown in the right panel of Fig|4] The temperature 
map reveals that subcluster A has a uniform temperature distri- 
bution in the range 4-5 keV with a drop in the inner 30" region 
to 3 keV, whereas subcluster B seems to have a temperature gra- 
dient, increasing from 2 keV in the surface brightness peak to 4 
keV. The abundances in every region chosen for spectral analy- 
sis show values in the range 0.3-0.5 Zq, with the exception of 
the central 0.5' region of subcluster A. In this region a single 
temperature fit returns a high metallicity of 0.76 ± 0.14 Zq. The 
regions coiTesponding to the surface brightness jumps described 
in Section [3] are consistent with being cold fronts, because in 
subcluster A the temperature increases from 4.2;°;^ keV in re- 
gion 2 to 5.0 ± 0.4 keV in region 3 and in the subcluster B from 
2.5 ± 0.2 keV in region 7 to SA^oi keV in region 8. 
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Fig. 4. Left panel: Regions used for spectral extraction to derive tlie temperature and abundance maps. The white circle corresponds to the ASCA 
aperture. Right panel: Black circles correspond to the values and la error bars for the temperature and abundances obtained with the IT model 
for each region shown in the left panel, discussed in the text in Section 1431 The dashed vertical lines indicate the regions adjacent to the surface 
brightness jumps of Figure|2l The red triangles correspond to the values for temperature and abundance obtained with a 2T model, as discussed in 
the text in Section l42l The higher temperature value at 5 keV has been fixed. 
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Fig. 5. Left panel: X-ray spectrum of the core of the main subcluster of Abell 2028 taken from a 0.5' aperture (corresponding to region 1 of 
Fig|4ll. Data from MOSl, MOS2, and pn are plotted in black, red, and green, respectively. The best-fit model obtained from a IT fit in the broad 
band 0.5-10 keV and residuals are also shown. 

Central Panel: MOS 1 spectrum from the inner core corresponding to region 1. Three different spectral models are overplotted: the IT model 
(green line) fitted in the broad 0.5-10 keV band; the IT model fitted in the 0.5-3 keV band (red line) and the IT model with the parameters 
obtained by the 0.5-10 keV fit, but with the abundance value kept fixed at the lower value obtained in the 0.5-3 keV band (blue line). 
Right panel: Same as in the left panel with the best-fit 2T model (with the higher temperature fixed at 5 keV) in the 0.5-10 keV band. 



4.2. The spectrum of the core of subcluster A 

We focus here on the spectrum of the core of the main subcluster, 
corresponding to region 1 in FigH] which is a circular region of 
0.5' radius, coiTesponding to 44 kpc. The single temperature fit 
performed in Section|4T|is a problematic fit because of residuals 
around the structure of the Fe-L shell, as can be seen in the left 
panel of Fig|5l and the high x^/dof = 243/168 highlights it 
further (see Table [T] for the results of all the fits discussed in 
the text for this region). The Fe-K line complex is not evident 
in the data, and a fit in an energy range (0.5-6.0 keV) which 
excludes the redshifted energy of the Fe-K line (6.2 keV) returns 
a metallicity of Z = 0.55 ± 0.10 Zq. If we restrict the energy 
band in the 0.5-3 keV range in an attempt to better model the 
Fe-L complex, we find a lower temperature (1.99 ± 0.09 keV 
with respect to 2.95 ± 0.13 keV) and abundance (Z = 0.32t°:E;^ 
Zq with respect to Z = 0.76 ± 0.14 Zq) compared to the broad 
band fit but with an improved fit, x^/dof = 178/148. Clearly 
this model underestimates the data at high energy when applied 



to the broad 0.5-10 keV band, as shown in the central panel of 
Fig|5]where we show the MOS 1 spectrum with the best-fit IT 
models in the broad and soft energy band overplotted with the 
shape of the Fe-L shell broad feature and residuals between the 
two fits in the inset. 

These results clearly indicate the need for a range of temper- 
atures for fitting the spectrum of the core of A2028 A and we 
therefore investigated fitting a two-temperature (2T) model to 
the data with the abundances tied between the two thermal com- 
ponents because the quality of the data does not allow relaxing 
this condition. The 2T fit does not constrain the temperature of 
the hotter component and the results of the fits obtained by fixing 
the temperature of the hotter component in the range 5-6 keV, as 
suggested by the temperature values in the regions of Fig|4l are 
listed in Table [T] Using a temperature of 6 ke V or 5 ke V for the 
hot component leads to statistically identical fits, with a temper- 
ature for the cooler component in the range 1.20-1.42 keV and 
metallicity in the range 0.53-0.75 Zq. The 2T model is a better 
description of the data with x^/dof = 204/167 corresponding 
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Fig. 6. X-ray spectrum of the inner 10" of the core of Abell 
2028 A (having subtracted the emission coming from an annulus 
with bounding radii 15"-30"). Data from MOSl, MOS2, and 
pn are plotted in black, red, and green, respectively. The best-fit 
model obtained from a IT fit and residuals are also shown. The 
pn spectrum taken from the 15"-30" annulus is shown in blue 
for comparison. 



4.3. The spectrum of the large ASCA aperture 

We extracted a spectrum from a c ircular region of 7' centered 
on RA 15:09:29 DEC H-07:34:20 (lHornerfl200l V which corre- 
sponds to the ASCA aperture (the white circle in Fig|4|i. We fit- 
ted a IT model and multi-temperature models as done for the 
core of A2028 A. As expected given the wide range of temper- 
atures revealed by the temperature map of Fig 01 the IT model 
fails to accurately model the spectra of this ASCA aperture which 
include the bulk of emission of the two subclusters. The multi- 
temperature models also do not fit the spectra adequately. But 
significantly we do not see the large variations in metallicity 
among different models as seen in the case of the core of Abell 
2028 A, and it is close to the average abundance in the various re- 
gions of this merging cluster as revealed by the two-dimensional 
map. There are enough counts to allow us to fit the spectrum 
just in the 2-10 keV energy band and to measure an abundance 
based just upon the Fe-K line. The results are consistent within 
the errors with all the previous determinations (0.35 ± 0.05 Zq 
obtained with a IT fit in the 2-10 keV band, to be compared with 
the upper value of the range, 0.43 ± 0.03 Zq obtained with the 
IT and GDEM models in the 0.5-10 keV band, and the lower 
value, 0.33 ± 0.04 Zq obtained with a 2T model in the 0.5-10 
keV band). 



to a significance of the additional component, calibrating the F- 
statistic using simulati ons of the null model (the IT model fit 
above) as suggested bv lProtassov et al.l (l2002h . of 99.99 %. 

We investigated other multi-temperature models to attempt 
to constrain the general differential emission measure (DEM) 
of the emission spectrum; in particular, following Buote et al. 
(I2003ah . we used a model in which the temperature distribu- 
tions is a Gaussian (GDEM) centered on the mean To with 
width (Tt (one additional parameter compared to the IT model) 
and a model where the temperature distribution is a power law 
(PDEM) with slope a and width of the temperature distribu- 
tion characterized by T^ax and T,nin (two additional parameters 
compared to the IT model). The abundances are tied between all 
the thermal components. The fits with these multi-temperature 
models are comparable to the one obtained with the 2T model, 
and they confirm the composite nature of the spectrum of the 
core. We also performed the fits using t he MEKAL plasma code 
dMewe et al.lll985l:lLieda"hr et al.|[T995 l) finding similar results. 

Given the evidence of a soft component with temperature in 
the range 1.3-1.5 keV and the surface brightness profile shown 
in Figl3] we made an attempt (bearing in mind all the caveats 
of extracting a region comparable in size with the half energy 
width of the PSF) to extract a spectrum from a circular region 
of 10" radius using as background a spectrum extracted from 
an annulus with bounding radii 15"-30", both centered on the 
emission peak. The purpose of this attempt is to investigate the 
excess emission that might be associated with the presence of 
a cool corona. The use of the local background very close to 
the excess emission should subtract the hot thermal emission of 
the surrounding ICM thereby revealing the cooler emission of 
the corona. The spectrum of the inner 10" is well fitted by a IT 
model with kT = 1.46 ±0.10 keV. We find the same result if we 
simultaneously fit the two spectra and we take into account the 
spreading of the cooler X-rays from the very center due to the 
finite PSF by means of th e cross-talk modifica tion of the ARE 
generation software (e.g. ISnowden e t al. 2008). The spectrum 
with the best-fit model is shown in Fig |6] together with the pn 
spectrum of the 15"-30" annulus to represent the hotter ICM in 
which the corona is embedded. 



5. Discussion 

5. 1. The inverse iron bias at work in Abell 2028 

The XMM-Newton observation of Abell 2028 clearly shows its 
complex structure by revealing two merging subclusters, a more 
luminous SW component (A2028 A) and a less luminous NE 
component (A2028 B). Associated with the brightest central 
galaxy (BCG) of A2028 A, GIN 416, there is indication of a cool 
corona as suggested by the ex cess surface brig htness profile and 
a cooler spectral component (ISun et al ] 120071) . The BCGs that 
able to retain their corona are usually associated with merging 
clusters, such as the first observed objects of this class, the BCGs 
in the Coma cluster ( Vikhlinin et al. 2001a). 

It is therefore clear that, within the beam of the ASCA re- 
gion us ed to derive the abu ndances presented for A2028 in the 
work of lfiaumgartner et al.l (12005). severa l components w i th dif- 
ferent temperatures are present. Work by Mazzotta et al. (2004|) 
and Vikhlinin (2006) pointed out how the temperature obtained 
by a fit with a single-temperature plasma model can be sig- 
nificantly different from commonly used definitions of average 
temperature like th e emission- weighted te mperature. However, 
as pointed out by iMazzotta et al.l (|2004), the spectrum of a 
multi-temperature model cannot be reproduced correctly by a 
single-temperature thermal model. On the other hand, a multi- 
temperature thermal emission may be well-fitted by a single- 
temperature model only because of the convolution with the in- 
strument response, Poisson noise and background noise. Even 
more difficult to gauge is the concept of the spectroscopic-like 
abundances, i.e. the measured abundance using a single tem- 
perature fit to a multi-temperature spectrum with possible dif- 
ferent abundances in different temperature components, because 
it depends strongly on the correct modeling of the temperature 
structure and on the limitations of current CCD detectors, which 
cannot spectrally separate emission from different components 
in a precise way. Besides the early discovery of the Fe bias 
(Buote 2000b), only recently work has been done to investigate 
this is sue (iMazzotta et al.ll2004tlBalestra et al]|2007l:lRasia et al.1 
I2008h . 
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Table 1. Results of the fits using single and multi-temperature models (see text for details of the models ) in different energy bands 
for the EPIC spectrum of the core (circle of 0.5' radius) of A 2028 A. 



Model/Band 


fcTo 


Z 


Normo 




Nornii 


CTT/width 




Notes 


(keV) 


(keV) 




(lO^^cm^^) 


(keV) 


(lO^^'cm^^) 


(keV) 






IT 0.5-10 


2.95 ±0.13 


0.76 ±0.14 


1.70 ± 0.09 








243/168 


Full spectrum 


IT 0.5-6 


2.59 ±0.10 


0.55 ±0.10 


1.93 ±0.09 








234/165 


no Fe-K 


IT 0.5-3 


1.99 ± 0.09 


0.321HI 


2.12 ±0.09 








178/148 


soft band Fe-L 


2T 0.5-10 
2T 0.5-10 
GDEM 0.5-10 


5.0 (fixed) 
4.0 (fixed) 

'^•''^-0.22 


0.53in? 
7ts+0-i3 

^- '>J-0.08 


1 42+0.06 
1.69 ±0.13 


i.48l«:?I 

1 97+0.07 
^•^ ' -0.03 


o.65l°:;| 

r, 97+0.02 


i.67l«:^^ 


204/167 
204/167 
202/167 




PDEM^ 0.5-10 


0.44 ± 0.08 


1 c;c:+0.45 


1.36 ±0.12 






5.8ltr7? 


0.83t!;:ig 197/166 





^ For PDEM Tq corresponds to Tmin and the width to Tmax — T, 




Fig. 7. Left panel: Exposure-corrected ASCA GIS image of Abell 2028. The regions inside the solid green 
lines (R=7') were extracted for the source spectrum analyzed in Baumgartner et al. (2005). Figure taken from 
|http://asd.gsfc.nasa.gov/Donald.Homer/acc/html/ABELL.2028.83040050.html 

Right panel: Comparison of the X-ray spectra of the large ASCA aperture from the ASCA SIS and 1 detectors (black and red, respectively) and 
from the EPIC pn (green). 



In the work by iRasia et al] ( l2008l) . who analyzed six 
simulated galaxy clu sters processed thro ugh an X-Ray Map 
Simulator (X-MAS, i GardinietaD |2004 that allows mock 
XMM-Newton MOS 1 and MOS 2 observations to be created, 
a systematic overestimate of iron for systems in the 2-3 keV 
range was found. By simulating a spectrum resulting from the 
combination of plasmas at 2 and 3 keV, they found fitted iron 
abundances higher with respect to the input value up to 40% in 
low S/N spectra (but still with 4500 counts). They also found the 
very interesting fact that the overestimate depends inversely on 
the number of counts, i.e. low S/N spectra have a greater bias. 
Projection and low resolution effects can create a complex tem- 
perature structure, averaging different temperatures in the same 
radial bin. Simionescu et al. (2009) give further support to this 
interpretation analyzing a 3' aperture in the center of the cool 
core cluster Hydra A finding that a GDEM model is a better fit 
to the spectrum and that it returns a best-fit Fe abundance lower 
than the one obtained from a fit with a single-temperature model 
of the full spectral band. (Lower values are derived when fitting, 
in particular, just the soft band.) 

The statist i cal qu ality of either the simulated spectra 
by iRasia et aP (l2008l) or the Hydra A spectra presented in 



ISimionescu et al. I (l2009l) are good, and in particular the Fe-K 
feature is clearly present. Many of the objects in the ASCA sam- 
ple did not have sufficient signal-to-noise spectra for the Fe- 
K line emission to be present with the same statistical signif- 
icance also because Baumgartn er et al.l (l2005h did not include 
very bright clusters to not unduly bias the results of their stack- 
ing analysis. This is indeed the case for the XMM-Newton spec- 
trum of the core of A 2028 A or the ASCA spectra of the 7' 
large aperture of the cluster (see right panel of Fig|7]l. With the 
aid of the simulations presented in appendix|A] we suggest that 
the inverse-Fe bias is mainly driven by the behavior of the the 
Fe-L complex strength which falls rapidly with increasing tem- 
perature above roughly 3 keV: in a multi-T spectrum where the 
average T is near 3 keV, but there are temperature components 
below and above 3 keV, most of the the Fe L emission of this 
multi-T spectrum comes from the lower-T components in the 
spectrum. However, when one fits a (wrong) single-T model, the 
single temperature will be sufficiently high that the Fe L lines 
will be weaker, and thus the model will have to compensate by 
increasing the Fe abundance above the true value. This is also 
shown by the central panel of Fig|5] That the Fe-K line is not 
necessary for the inverse iron bias to work can also be demon- 
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Fig. 8. Left panel: SDSS color image of the field of Abell 2028 Right panel: The RASS field of Abell 2028 highlighting in magenta the other 
members of the supercluster as Abell 2033 and Abell 2029. 



strated by simulating IT spectrum with the characteristics of the 
fit in the 0.5-3 keV band for the spectrum shown in Section I4r2l 
i.e. a single thermal spectrum with kT — 2.0 keV and Z = 0.38 
Zq . If we force the temperature in the model fit to be fixed at 3 
keV (as in the case of a broad band fit where higher temperature 
components are present in the spectrum), the abundance will be 
biased high at a value of Z = 0.68 Zq. The thermal structure 
of a merging cluster is different from the very central regions of 
a cool core cluster and the main driver of the inverse Fe-bias in 
A2028 A is the cool component of the corona, which also jus- 
tifies the simulations of Appendix lAl and fits with a 2T model, 
which seems a good physical description of a corona embedded 
in the hot intra-cluster gas. More complicated multi-temperature 
models, like the GDEM model, may provide on the contrary a 
more physical interpretation of the spectrum of the central re- 
gions of cool core clusters, as suggested bv ISimionescu et aP 
d2009h . 

The inverse-iron bias is, however, not present in the XMM- 
Newton EPIC spectra of the large 7' aperture. The S/N of this 
spectrum is higher and the presence of a cool component is 
less imp ortant. This reinforces the fact already pointed out by 
iRasia et al. (2008) that the bias depends inversely on the number 
of counts present in the spectrum. It is beyond the scope of this 
paper to investigate the reasons for the anomalous abundances 
directly in the ASCA spectra but the presence of the inverse iron- 
bias in those low S/N spectra (obtained with a lower effective 
area instrument) is a reasonable suggestion. The XMM-Newton 
results for A2028 show that the abundances in this cluster are 
not anomalously high. 

The presence of the direct u"on bias (lBuotdl2000bl) and of 
the inverse iron bias described in this paper show the danger of 
interpreting abundances measured with single-temperature fits 
to spectra extracted from large physical apertures (and of rela- 
tively low S/N), where there might be a great deal of temperature 
structure. This is particularly true for the combination of current 
telescopes and CCD detectors, which give a lot of weight to the 
soft energies and therefore to the shape of the Fe-L complex. 
Caution should be paid, for example, when interpreting global 
abundances measured in clusters out to a redshift of z ~0.5 be- 
cause the Fe-L shell of low-temperature components (down to 



^1 keV) is within fitting bandpasses with lower energy ends of 
0.5-0.6 keV. 

5.2. The dynamical state of Abell 2028 

The X-ray morphology of the two subclusters associated with 
Abell 2028, with their tails and cold fronts, clearly sug- 
gest s an off-axis merger as seen in various simulations (e.g. 
Ric ker & Sarazini r2001). It is therefore interesting to look for 
other information about the dynamical state of the cluster at dif- 
ferent wavelengths. Abell 2028 in the optical band is an elon- 
gate d structure as shown by galaxy iso-density images (see Fig. 
8 of iTrevese et al.|[T997i) or by wavelet analysis (see the panels 
on sc ales of 750 kpc and 527 kpc of Fig. 13 in Flin & K rvwultl 
'2006"). On a scale of 258 kpc, iFlin & Krvwultl (|2006 ) detected 
substructures in the 2D galaxy distribution. These structures 
are coincident with the X-ray subclusters. The surface bright- 
ness peak of A 2028 A is coincident with the galaxy GIN 416, 
the BCG of the subcluster in the optical. The peak of A 2028 
B is coincident with the pair of galaxies LEDA 140479 and 
2MASX J15092098H-0738225, which seem to be interacting (see 
the left panel of Fig|8]l. The velocities of these three galaxies are 
cz = 23148 km s^^, cz = 23011 km s-\ and cz = 23135 
km s~^, suggesting a merger very close to the plane of the 
sky. The definitive confirmation of the presence of a corona has 
to wait for Chandra data with its better PSF; however further 
support is given by GIN 416 also being a radio source, NVSS 
150928H-073324, with = 4.1 x 10^3 W Hz^ (neglecting the 
small K-correction term), which is consistent with the idea that 
a radio active BCG should at least have a corona to be in such a 
state, because there is no large-scale cool core (Sun 2009) . 

The galaxy distribution is elongated in the direction of the 
closest neighbor Abell 2033. Abell 2028 and Abell 2033, to- 
gether with Abell 2029 and Abell 2066, for m the s uper-cluster 
number 154 in the catalogue of lEinasto et al.l (12001 1) . It is likely 
that another cluster can be added to the membership of the su- 
percluster: the brig ht /JO^Ar source IRXS J150935.9H-084605 
(IVoges et al.lll999|) is associated with a clustering of 25 SDSS 
galaxies, with measured spectroscopic redshifts, at the redshift 
of the supercluster (a Gaussian fit to the histogram of the galax- 
ies with measured redshift gives cz = 23948 km s^^ corre- 
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Fig. 9. Left panel: Summed MOS and pn surface brightness profile across the cold front in the Al sector of subcluster A shown in Fig|2l together 
with the best-fit model of a power law to describe the outer ICM and a /3-model to describe the density inside the edge. Right panel: Same as the 
left panel for the cold front in sector B 1 of subcluster B. 



spending to a redshift of z=0.0799, with a velocity dispersion 
cr = 574 ± 198 km s^i). The brightest galaxy in the field, SDSS 
J150936.33H-084632.8, is a red elliptical whose photometric red- 
shift measured by template fitting is 0.0788 ± 0.009. 



The subclusters of A 2028 are therefore merging to produce 
a more massive cluster belonging to this large-scale superclus- 
ter. If we coarsely indicate the direction of merging as indicated 
by the "tails" of the two subclu sters (and by the m erging cold 
fronts perpendicular to them, see iMarkevitch & Vikhlinin,2007. 
and references therein), and the direction of the super-cluster as 
indicated by the four members in the right panel Fig|8] A2028 
A and B are merging at ^ 60° direction with respect to the axis 
of the super-cluster In order to bound in a quantitative way the 
motion of the two subclusters we determine the jump in gas den- 
sity at the fronts, assuming that the gas density profile is de- 

-3/20 



scribed by the standard /3 model n = no 1 + {r/vcY 
inside the cold front (r > rcj where Vcf is the cold front ra- 
dius) and by a power law n = rii (r)~^ outside the cold front 
ir > Tcf). We then calculate the corresponding projected sur- 
face brightness. We fix the cold front radius at Vcf = 65" for 
the subclump A and Tc/ = 65" for subclump B. We fit the outer 
part of the surface brightness profile (summing the MOS and 
pn profiles shown in Fig|2|i to set the external component pa- 
rameters, and we successively derive the best fit parameters of 
the innermost part. In Fig|9]we plot the surface brightness pro- 
files with the best-fit models. We find a density discontinuity 
at the cold front in subcluster A of Uin/nout — 3.11 ± 0.57, 
which combined with the temperature values in the regions in- 
side and outside the front (see Sectior l^TI ). leads to a pressure 
ratio Pin/ Pout = 2.61 ± 0.53. In subclump B, we find a den- 
sity discontinuity at the cold front nin/riout = 5.86 ± 1.17, 
which combined with the temperature values, l eads to a pres- 
sure ratio Pin/ Pout = 4.31 ± 0.98. Following fVikhlinin et aP 
(,2001b). these pressure ratios correspond to a Mach number 
M = 1.19to:i8 for subclump A and X = 1.61^1^:22 for sub- 
dump B respectively. 



6. Conclusions 

We analyzed an XMM-Newton observation of Abell 2028, a ~ 4 
keV cluster previously known with ASCA data to display anoma- 
lously high global abundances. The main results of our work can 
be summarized as follows: 

- Abell 2028 is actually composed of two subclusters in the 
process of merging, along the axis of a supercluster. The 
evidence suggesting this interpretation is the cometary-like 
structure of the two subclusters with cold fronts at the lead- 
ing edges of the two structures. The brightest galaxy of the 
main subcluster is likely hosting a cool corona; 

- The abundance of the clusters, as derived from a tempera- 
ture and abundance map, is consistent with the local abun- 
dance of more massive clusters. The abundances in every re- 
gion chosen for spectral analysis show, with the appropriate 
temperature modeling, values in the range 0.3-0.5 Zq, much 
lower than the previous ASCA global me asurement of 0.77 
Zq (in units of iAnders & Grevesselll989l) : 

- We suggest that the inverse Fe-bias, i.e. the overestimate of 
iron abundance when fitting multi-temperature spectra with 
average spectroscopic temperature around 3 keV with a sin- 
gle temperature model, is likely the source of the high abun- 
dances in 3-4 keV local clusters as already discussed in the 
literature. We suggest that the behavior of the Fe-L shell lines 
explains the bias. The Fe-L complex strength falls rapidly 
with increasing temperature above roughly 3 keV: in a multi- 
T spectrum where the average T is near 3 keV, but there are 
temperature components below and above 3 keV, most of 
the the Fe L emission of this multi-T spectrum comes from 
the lower-T components in the spectrum. The best-fit single 
temperature will be high enough that the Fe L lines will be 
weaker, so the model will have to compensate by increasing 
the Fe abundance above the true value. 

More definitive results will be obtained with our ongoing 
project of observing local clusters in the 3-4 keV range. More 
quantitative results for the dynamical state of A 2028 and for the 
cold fronts and corona detected in the XMM-Newton data will be 
possible with a dedicated Chandra follow-up. 
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Appendix A: Analysis of simuiated X-ray spectra 

In this section we investigate the inverse Fe bias when fitting 
simulated multi-temperature models with a single-temperature 
plasma model. Following a well-established approach we as- 
sumed a double temperature model (2T) to si mulate a plasma 
with substantial temperature structure (e.g. iBuotd l2000bt 
Balestra et al. 2007; Rasia et al. 2008). Even though it may be 
argued that a model with a more continuous distribution of emis- 
sion measure could be more physically motiva ted (see for exam- 
ple th e discussion for the core of Hydra A in ISimionescu et al.l 
I2009 ?). it is difficult to disctinguish it with respect to a 2T model; 
this also seems to be realistic on physical grounds given the pos- 
sible interpretation of a cool corona embedded in an ambient 
hotter plasma we put forward to explain the core of Abell 2028 
A, which we used as a reference starting point for our simula- 
tions. We perform several simulations of spectra with different 
assumptions, as described in detail below, and explore the possi- 
ble conditions that can potentially affect the distribution of best- 
fit values of kT and, most important, of the metallicity Z. 

We performed a set of simulations of clusters with two tem- 
perature components (specified in Table lA.lb in the range rel- 
evant for the core of Abell 2028 A, i.e. 1.3-2 keV for the cold 
component and 4-6 keV for the hot component with an emission 
measure ratio (the normalization of the APEC model) ranging 
from 0.1 to 1 .0 with steps of 0. 1 between the two temperature 
components. The starting dominant component is listed as first, 
with an initial normalization (in XSPEC units) of 1.5 x 10^''. 
The total number of counts change as a function of the ratio 
EMcoid/EMhot and ranges from ^1100 counts and ^2600 for 
MOS and pn in the 0.5-10 keV band in the simulations with 
EMcoid/EMhot = 0.1 to -1800 and -4600 for MOS and pn 
in the simulations with EMcoid/EMhot = 1-0. For each set of 
simulations, we performed 20 realizations and listed the results 
of the fits with a IT model in the various table of results as me- 
dian and mean absolute deviation from the median. We simu- 
lated spectra with the Galactic column density and redshift of 
A 2028 with an exposure time of 50 ks. We also investigated 
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the distribution of best-fit abundances in the presence of a gradi- 
ent in the abundance profile, with a higher abundance associated 
with the lower temperature component. We quote and plot the 
results separately for the three EPIC detectors to investigate the 
presence of effects arising by the difference in effective area. 

It can be seen by the plots of the results for the various simu- 
lations (from Fig. lA.ll to lA.6b . which are also listed in Table lA.2l 
that, as the contribution of the low temperature component in- 
creases, the best-fit temperature of the single component model 
shifts to lower values. The distribution of the best fit abundance, 
either in cases with the same input abundance for the two compo- 
nents or in cases with an higher abundance in the lower temper- 
ature component, is always biased to higher values that a naive 
emission measure weighting would predict, up to ~ 50% of the 
median value compared to the input value in some cases. Even 
more striking is the large scatter in the distributions in particu- 
lar for low values of the ratio EMcoid/EMhot and for the MOS 
detectors, again pointing to the relevance of the Fe-L line in this 
bias: the fitting of the feature is key for the abundance deter- 
mination, and Poisson noise is limiting for lower effective area 
instruments. Rather extreme values may occur: we show from 
example a MOS 1 si simulation when EMcoid/EMhot = 0.3 
(see Fig |A.7b : the best-fit abundance is 0.72 ± 0.23 Zq with a 
statistically good fit (x^/dof < 1) with just a single temperature 
(we recall that the input abundance value for both components 
is 0.3 Zq). It is instructive to plot this example together with 
another realization of the same model, with a best-fit abundance 
of O.OT^Q gy Zq, because it shows the importance of the fitting 
of the spectral region of the Fe-L lines (see Fig lA.Sb . In both 
cases the best fit temperatures are around ^ 3 keV. That the role 
played by the contribution of the Fe-K line is not relevant in the 
inverse Fe bias is also suggested by the simulation s6, where we 
only use a bremsstrahlung (i.e. the abundance was set to zero) 
for the hot component. Already for EMcoid/EMhot = 0.4 the 
median best fit abundance is equal to the input abundance of the 
cold component. 

The trends shown here do not change if one instead selects 
spectra with the same number counts as the ratio EMcoid/EMhot 
increases. As an example we show the results in Fig |A.9| of the 
simulation si for MOS2 for three different number counts for 
all the EMcoid/EMhot ratios, i.e. 1 x 10^, 5 x 10^ and 1 x lO^* 
counts in the 0.5-10 keV band. 

Another result coming from the investigation of the simu- 
lations is that restricting the fitting energy band to just fit the 
shape of the Fe-L bump (for example in the 0.5-3 keV energy 
band) is no guarantee of recovering the abundance of the low- 
temperature component: the fitted continuum is too high because 
of the contribution of the high-temperature component and the 
abundance thus fitted is too low (see the green and blue pn spec- 
tra of Fig|6]for a visual representation). 



Table A.l. Input values for the spectral simulations of two tem- 
perature components as described in the Appendix 



Bim" 


kT ^ 


7 


si 


4+1.3 


0.3 


s2 


4+1.3 


0.2+0.4 


s3 


4+2 


0.2 


s4 


4+2 


0.2+0.4 


s5 


5+2 


0.2+0.4 


s6 


6+2 


0.0+0.4 



Notes: ^ simulations identification number; input temperature of the 
hot and cold component; input iron abundance; XSPEC fixed 
norm of the first temperature component (cm"""); 




0.2 0.4 0.6 0.8 1 



Fig. A.l. Median of the best-fit temperature and abundance dis- 
tribution as a function of the ratio of the cold and hot emission 
measure, for the simulation si. The green dots refer to Ml re- 
sults, the red dots to M2 and the blue dots to pn. The dashed 
lines refer to the input values of temperature and abundance for 
the two components. The error bars correspond to the mean ab- 
solute deviation. 





■ * * i 

* * ' » 















0.2 0.4 0.6 0.8 1 



Fig.A.2. Same as in Fig lA.ll for the simulations s2. Note the 
two lines in the abundance panel which refer to the two different 
input abundances in the two components: the higher one (0.4 
Zq) in the cold component and the lower one (0.2 Zq) in the 
hot component. 
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Fig. A.3. Same as in Fig lA. ll for the simulations s3. 
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Fig. A.4. Same as in Fig lA. ll for the simulations s4 
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Fig. A.5. Same as in Fig lA. ll for the simulations s5 
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Fig. A.6. Same as in Fig lA. II for the simulations s6 



1^^ 





Energy (keV) 

Fig. A.7. The MOS 1 spectrum and best fit IT model with resid- 
uals for a realization of the simulation si with EMcoid/EMhot = 
0.3. The input model consists of two plasmas with kTcoid = 1-3 
keV and kT^ot — 4 keV with an abundance of 0.3 Zq in both 
components, the single-temperature model best fit values are 
kT = 3.02 ± 0.22 keV and Z = 0.72 ± 0.23 Zq. 

Comparison M! s2 EM^,^,|/EM^_,:^0.3 



ft 



Fig.A.8. Comparison of the MOS 1 spectra for two different 
realizations of the simulation si with EMcoid/EMhot — 0.3: the 



black dots refer to the realization of Fig lA. 7l and the red dots to 
a realization with best fit abundance of 0.07 Zq. 
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Fig. A.9. Median of the best-fit temperature and abundance dis- 
tribution as a function of the ratio of the cold and hot emission 
measure, for the simulation si for the detector MOS2, with the 
same number counts for all the EMcoid/EMhot ratios. The red 
dots refer to spectra with 1 x 10'^ counts in the 0.5-10.0 keV 
band, the green ones to spectra with 5 x 10"^ counts, and the blue 
to spectra with 1 x 10^ counts. The dashed lines refer to the input 
values of temperature and abundance for the two components. 
The error bars correspond to the mean absolute deviation. 
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Table A.2. Fitted values with a IT model in the 0.5-10 keV band for the three EPIC detectors of the simulations listed in Table I A. fl 





Ratio 0. 1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 












.si 


Ml 










kT 


3.34 ± 0.22 


3.11 ± 0.16 


2.92 ± 0.17 


2.74 ± 0.14 


2.59 ± 0.14 


2.48 ± 0.14 


2.39 ± 0.15 


2.09 ± 0.17 


1.99 ± 0.13 


1.98 ± 0.12 


Z 


0.26 ± 0.14 


0.43 ± 0.12 


0.40 ± 0.14 


0.41 ± 0.15 


0.35 ± 0.11 


0.35 ± 0.08 


0.36 ± 0.07 


0.36 ± 0.10 


0.31 ± 0.06 


0.35 ± 0.07 












si 


M2 










kT 


3.42 ± 0.29 


3.10 ± 0.27 


2.88 ± 0.20 


2.74 ± 0.17 


2.53 ± 0.08 


2.44 ± 0.11 


2.32 ± 0.15 


2.07 ± 0.15 


2.05 ± 0.10 


1.97 ± 0.07 


Z 


0.22 ± 0.23 


0.31 ± 0.15 


0.34 ± 0.13 


0.35 ± 0.08 


0.31 ± 0.07 


0.36 ± 0.08 


0.37 ± 0.07 


0.37 ± 0.06 


0.35 ± 0.07 


0.31 ± 0.06 












si 


pn 










kT 


3.47 ± 0.16 


3.15 ± 0.16 


2.84 ± 0.14 


2.59 ± 0.13 


2.47 ± 0.09 


2.37 ± 0.12 


2.30 ± 0.14 


2.08 ± 0.13 


2.01 ± 0.07 


1.96 ± 0.06 


Z 


0.32 ± 0.06 


0.30 ± 0.07 


0.37 ± 0.06 


0.42 ± 0.07 


0.35 ± 0.06 


0.35 ± 0.06 


0.35 ± 0.06 


0.37 ± 0.05 


0.36 ± 0.05 


0.35 ± 0.04 












s2 


Ml 










kT 


3.51 ± 0.21 


3.11 ± 0.24 


2.84 ± 0.22 


2.61 ± 0.17 


2.43 ± 0.11 


2.07 ± 0.14 


2.02 ± 0.09 


1.97 ± 0.06 


1.90 ± 0.09 


1.86 ± 0.13 


Z 


0.24 ± 0.13 


0.33 ± 0.17 


0.39 ± 0.11 


0.40 ± 0.13 


0.44 ± 0.08 


0.35 ± 0.07 


0.35 ± 0.07 


0.33 ± 0.07 


0.37 ± 0.04 


0.33 ± 0.05 












s2 


M2 










kT 


3.35 ± 0.24 


3.00 ± 0.17 


2.67 ± 0.22 


2.50 ± 0.13 


2.44 ± 0.17 


2.24 ± 0.18 


2.01 ± 0.10 


1.98 ± 0.08 


1.93 ± 0.08 


1.89 ± 0.10 


Z 


0.32 ± 0.12 


0.31 ± 0.09 


0.32 ± 0.11 


0.37 ± 0.08 


0.47 ± 0.12 


0.36 ± 0.08 


0.35 ± 0.07 


0.34 ± 0.06 


0.34 ± 0.07 


0.37 ± 0.07 












s2 


pn 










kT 


3.36 ± 0.18 


3.08 ± 0.16 


2.65 ± 0.15 


2.49 ± 0.08 


2.40 ± 0.10 


2.28 ± 0.17 


2.02 ± 0.14 


1.97 ± 0.06 


1.89 ± 0.03 


1.88 ± 0.03 


Z 


0.34 ± 0.11 


0.36 ± 0.05 


0.36 ± 0.04 


0.36 ± 0.04 


0.39 ± 0.05 


0.36 ± 0.05 


0.38 ± 0.05 


0.36 ± 0.04 


0.37 ± 0.04 


0.37 ± 0.04 












s3 


Ml 










kT 


3.71 ± 0.30 


3.43 ± 0.29 


3.25 ± 0.26 


3.15 ± 0.23 


3.16 ± 0.22 


3.19 ± 0.23 


3.09 ± 0.17 


3.01 ± 0.15 


2.88 ± 0.16 


2.79 ± 0.15 


Z 


0.28 ± 0.14 


0.25 ± 0.10 


0.22 ± 0.09 


0.24 ± 0.09 


0.22 ± 0.12 


0.26 ± 0.07 


0.26 ± 0.12 


0.26 ± 0.10 


0.26 ± 0.08 


0.24 ± 0.07 












s3 


M2 










kT 


3.55 lb 0.28 


3.44 ± 0.21 


3.40 ± 0.21 


3.33 ± 0.22 


3.24 ± 0.20 


3.16 ± 0.18 


3.12 ± 0.16 


3.07 ± 0.13 


2.93 ± 0.12 


2.83 ± 0.15 


Z 


0.23 ± 0.11 


0.30 ± 0.15 


0.29 ± 0.14 


0.27 ± 0.11 


0.28 ± 0.09 


0.32 ± 0.13 


0.25 ± 0.09 


0.25 ± 0.09 


0.30 ± 0.11 


0.22 ± 0.11 












s3 


pn 










kT 


3 58 + 18 


3 39 + 


3 93 + 16 


3.24 ± 0.11 


3 16 + 011 


3 05 + OQ 


3 01 +01? 


2.91 ± 0.12 


? 9? + 1 3 


9 70 + n 08 


Z 


0.21 ± 0.09 


0.29 ± 0.09 


0.25 ± 0.08 


0.24 ± 0.10 


0.27 ± 0.06 


0.25 ± 0.06 


0.21 ± 0.04 


0.20 ± 0.05 


0.24 ± 0.05 


0.25 ± 0.05 












s4 


Ml 










kT 






tj.^fj _i_ W.-Lt: 


3 07 + 1 8 


3 05 + ?1 


2.94 ± 0.14 


9 QQ _|_ n 1 f; 


2.84 ± 0.11 


2.82 ± 0.14 


9 71 -U n 1 


Z 


0.19 ± 0.16 


0.34 ± 0.14 


0.33 ± 0.15 


0.31 ± 0.15 


0.32 ± 0.08 


0.37 ± 0.10 


0.41 ± 0.11 


0.42 ± 0.09 


0.45 ± 0.08 


0.41 ± 0.10 












s4 


M2 










kT 


O . U -L _i_ V/ . ^ O 






1 4 _|_ n on 


3 15 + 016 


3 07 + 1 7 


2.94 ± 0.14 


? 86 + 1 3 




9 74 + n 1 5 


Z 


0.25 ± 0.18 


0.33 ± 0.19 


0.35 ± 0.17 


0.33 ± 0.14 


0.34 ± 0.09 


0.40 ± 0.07 


0.45 ± 0.14 


0.38 ± 0.12 


0.40 ± 0.10 


0.39 ± 0.08 












s4 


pn 










kT 


3.54 ± 0.22 


3.32 ± 0.11 


3.20 lb 0.12 


3.16 lb 0.14 


3.09 ± 0.16 


2.85 ± 0.10 


2.90 ± 0.07 


2.81 ± 0.09 


2.72 ± 0.10 


2.63 ± 0.08 


Z 


0.28 ± 0.09 


0.32 ± 0.08 


0.38 ± 0.07 


0.38 ± 0.07 


0.34 ± 0.08 


0.35 ± 0.06 


0.38 ± 0.06 


0.37 ± 0.06 


0.41 ± 0.07 


0.40 ± 0.05 












s5 


Ml 










kT 


4.61 ± 0.47 


3.92 ± 0.28 


3.88 ± 0.27 


3.59 ± 0.20 


3.40 ± 0.18 


3.22 ± 0.17 


3.27 ± 0.15 


3.11 ± 0.14 


3.02 ± 0.15 


3.00 ± 0.13 


Z 


0.26 ± 0.21 


0.36 ± 0.18 


0.32 ± 0.18 


0.37 ± 0.16 


0.44 ± 0.11 


0.44 ± 0.11 


0.44 ± 0.12 


0.44 ± 0.10 


0.46 ± 0.10 


0.43 ± 0.12 












s5 


M2 










kT 


4.53 ± 0.38 


4.15 ± 0.25 


3.83 ± 0.19 


3.60 ± 0.21 


3.36 ± 0.25 


3.32 ± 0.23 


3.32 ± 0.22 


3.17 ± 0.17 


3.15 ± 0.15 


3.09 ± 0.16 


Z 


0.25 ± 0.18 


0.42 ± 0.11 


0.38 ± 0.14 


0.42 ± 0.18 


0.44 ± 0.11 


0.49 ± 0.13 


0.47 ± 0.14 


0.41 ± 0.10 


0.47 ± 0.12 


0.48 ± 0.09 












s5 


pn 










kT 


4.41 ± 0.31 


4.02 ± 0.20 


3.82 ± 0.18 


3.63 ± 0.20 


3.48 ± 0.16 


3.34 ± 0.14 


3.18 ± 0.11 


3.13 ± 0.13 


3.04 ± 0.05 


2.93 ± 0.09 


Z 


0.27 ± 0.13 


0.36 ± 0.08 


0.37 ± 0.08 


0.32 ± 0.08 


0.43 ± 0.08 


0.42 ± 0.07 


0.41 ± 0.09 


0.42 ± 0.08 


0.47 ± 0.08 


0.43 ± 0.05 












s6 


Ml 










kT 


5.23 ± 0.48 


4.87 ± 0.42 


4.26 ± 0.29 


4.05 ± 0.24 


3.78 ± 0.29 


3.36 ± 0.23 


3.30 ± 0.18 


3.20 ± 0.18 


3.12 ± 0.18 


3.12 ± 0.16 


Z 


0.10 ± 0.18 


0.24 ± 0.15 


0.33 ± 0.18 


0.39 ± 0.17 


0.33 ± 0.13 


0.37 ± 0.11 


0.41 ± 0.16 


0.40 ± 0.09 


0.40 ± 0.11 


0.39 ± 0.11 












s6 


M2 










kT 


5.09 ± 0.47 


4.47 ± 0.28 


4.13 ± 0.30 


3.91 ± 0.27 


3.74 ± 0.22 


3.64 ± 0.30 


3.30 ± 0.23 


3.16 ± 0.18 


3.14 ± 0.20 


3.02 ± 0.16 


Z 


0.13 ± 0.16 


0.27 ± 0.18 


0.36 ± 0.21 


0.34 ± 0.13 


0.36 ± 0.13 


0.38 ± 0.10 


0.40 ± 0.09 


0.39 ± 0.09 


0.44 ± 0.11 


0.47 ± 0.10 












s6 


pn 










kT 


5.30 ± 0.34 


4.92 ± 0.33 


4.51 ± 0.25 


4.11 ± 0.21 


3.83 ± 0.26 


3.57 ± 0.17 


3.33 ± 0.15 


3.14 ± 0.09 


3.02 ± 0.13 


3.11 ± 0.14 


Z 


0.05 ± 0.07 


0.13 ± 0.10 


0.18 ± 0.12 


0.24 ± 0.06 


0.24 ± 0.09 


0.35 ± 0.07 


0.36 ± 0.07 


0.37 ± 0.08 


0.35 ± 0.05 


0.38 ± 0.06 



